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A  METHOD  FOR  ESTIMATING  DANGER  AREAS  DUE 
TO  RICOCHETING  PROJECTIIES 


ABSTRACT 

Tne  rlc-chet  studies  of  Birkhoff  and  Hitchcock  are  summarized.  A  theory  is 
developed  for  ricochet  of  projectiles  off  inclined  plares.  and  this  theory  is 
applied  to  obtain  a  procedure  for  the  estimation  of  danger  areas  dt'“  t  -  ricochei.in£ 
projectii.es . 


I.  INTRODUCTION 


Ir.  several  instances  over  the  past  few  years  it  itsz  been  necessary  to  estimate 
the  danger  area  resulting  from  the  ricochet  01  projectiles  off  various  surfaces 
during  outdoor  training  exercises  or  testing  operations.  Of  particular  interest 
was  the  lateral  extent  of  the  danger  cone.  To  supply  these  estimates,  an  analyti¬ 
cal  nrooedure  was  worked  out  and  is  described  in  the  present  report. 

Based  upon  a  large  collection  of  experimental  data  on  ricochet  off  .land 
surfaces,  The  follow  -ng  empirica'1  formulas  were  published  by  G.  Birkhofi  in  19 4 j 
(BRL  Report  No.  535): 
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Where  (see  Figure  l):  9  =  arigle  of  impact,  |u|  -  magnitude  of  velocity  of 

impact,  0  -  angle  of  ricochet,  |v|  =  magnitude  of  velocity  of  ricochet,  and 

©c  =  the  critical  angle  for  the  projectile/ soil  (or  intact  material)  combination 
and  is  defined  as  that  angle  of  impact  above  which  ricochet  occurs  less  than  half 
the  time. 
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In  lyVf,  H.  P.  Hitchcock  published  a  large  collection  of  additional  experi¬ 
mental  data  on  ricochet  off  land  surfaces  (BRL  Report  Ilo.  (icy).  Considering  these 
additional  data,  D.  o .  Dunn  obtained  the  following  modified  Birkhoff  formulas : 
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3.5  -  2.s  -f-  (e  <  8c) 


(3) 

w 


Graphs  of  Equations  (5)  and  (4)  are  presented  in  Figure  2.  These  luodified 
Birkhoff  formulas  are  currently  used  at  BRL  for  computation  of  velocity  and  angle 
of  ricochet  from  velocity  and.  angle  of  impact.  The  critical  angle  ©c  is  a  function 
of  the  nose  shape,  impact  velocity,  hardness  of  the  projectile  nose  and  hardness 
of  the  ricochet  surface.  The  value  of  y,  is  obtained  by  searching  through  the 
ricochet  literature  for  an  experimental  value  corresponding  to  the  appropriate 
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sat  of  parameters.  Notice  in  Equations  (3)  ann  (h),  or  in  Figure  2,  that  both 

jv|  and  0  increase  with  ail  increase  in  6  —  assuming  given  values  of  the  impact 

C 

conditions  |uj  and  0.  Consequently,  for  computation  of  danger  zones,  9^  should 
he  chosen  as  large  as  realism  permits  in  order  to  avoid  underestimating  the  size 
of  the  danger  zone. 

XI.  THEORY  OF  RICOCHET  OFT  INCLINED  PLANES 

Consider  the  impact  of  a  given  projectile  on  an  inclined  plane  composed  of  o 
given  material.  We  assume  that  the  critical  angle,  9c,  for  the  projectile/ impe ■"+ 
material  combination  has  been  determined  (or  appropriately  chosen).  The  point  of 
impact  is  taken  as  the  origin  (0,  0,  0)  of  an  xyz  rectangular  coordinate  system, 
Figure  3*  The  y  axis  is  an  extension  of  the  radius-vector  from  the  center  of  the 
earth  to  the  origin  (0,  0,  0),  and  the  xz  plane  is  perpendicular  to  the  y  axis. 

-X 

The  x  axis  is  chosen  so  that  the  impact  velocity,  u,  of  the  projectile  is  contained 
in  the  xy  plane;  and  the  z  axis  is  then  chosen  so  that  the  xyz  system  will  be  right- 
handed.  Unit  vectors  along  the  x,  y,  and  z  axes  are.  denoted  fcy  i.  k,  respectively. 
Hence  we  have 


v +  V 


lul 


V  u2  +  u2 
X  y 


(5) 


(C) 


We  assume  that  u  and  u  ar<  known  from  the  output  of  a  two-dimensional  trajectory 
y  _  v 

confutation  terminating  at  (0,  0,  0).  Hence,  |u|  is  also  known.  The  orienta¬ 
tion  of  the  inclined  plane  is  specified  by  a  unit  vector  n  perpendicular  to  the 
inclined  plane  at  (0,  0,  0).  We  hare 


u 


f  rr\ 
\  1  / 


I  nl  =  -I  n2  +  r.2  +  n2  =  1 
xyz 


(S) 


The  two  components  nx  and  ny  will  be  arbitrarily  specified,  and  will  then  be 
determined,  except  for  sign,  from  Eqiation  (3).  Adopting  the  convention  that 
n  ^  0  (which  amounts  to  considering  half  of  a  symmetrical  distribution  cf  possible 
orientations  of  the  inclined  plans  --  i.e.  symmetrical  about  the  xy  plane)  we  have 
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The  orientation  of  the  inclined  plane  can,  of -course,  he  varied  by  varying 
and  n  .  For  any  particular  orientation,  our  known  quantities  are  u  ,  u  (and 

v  ^  ^  J 

hence  I  u  1  ),  n  ,  n  ,  n,  ,  and  6  .  In  terms  of  tnese  known  quantities,  we  must 

^  y  ^  v 

cal- -1 axe  (l)  the  magnitude,  |  v I  ,  of  the  velocity  v  of  ricochet,  (2)  the 
angle  a  that  the  vecxor  v  makes  with  the  xz  pine,  and  (3)  the  angle  -4.  that  the 
projection  of  v*  on  the  xz  plane  makes  with  the  x  axis  (see  Figure  3).  'The 
initial  conditions  (i.e.  initial  velocity  and  angle  of  elevation)  for  the  ricochet 
trajectory  are  |  v'  |  and  a,  ana  the  range  of  the  ricochet  trajectory  must  be 
plotted  so  as  to  form  the  angle  y  witn  the  x  axis. 

We  will  begin  by  determining  the  angle  of  impact,  0,  that  the  vector  u  makes 
with  the  inclined  plane  (see  Figure  3).  We  have  (see  Figure  4): 

u  •  n  =  lu|  \t\  cos(90°  +  ©) 

=  |u  i  (  -  sin  0)  (iO) 

_S  — 1  -A  1  u  u  ^ 

But,  since  u  =  u  i  +  u  j  and  since  n  =  ni+njxnk  ,we  have 
x  y  x  y  z 

u  .  n  =  u  n  +  u  n  (11) 

x  X  ^  V 

Equations  (lO)  and  (ll)  yield  the  following  solution  for  0  in  terms  of  the  known 
quantises  u^,  uy,  nx,  ny,  !~u|  : 


9  =  sin 


,  f  rvt  n  +  u  n  I 

L  lui  ) 


In  general  we  ha to  u  >  o.  n  <-  0.  u  <0.n  >0,  so  that  0  will  generally  turn 
x  x'y*y 

out  to  be  a  positive  angle. 

We  now  know  |uj  ,9,  and  9^  ;  and  so  we  can  compute  |v(  and  0  from  the 
modified  Birkhc»f  formulae,  Equations  (3)  ana  (4).  Of  course,  this  0  is  the  angle 

— k  -4 

that,  the  vector  v'  makes  with  the  inclined  plane.  Hence,  we  can  now  add  |vj 
and  0  to  our  list  of  known  quantities.  To  determine  the  angles  a  ana  4,  we  need 
to  find  the  components  of  the  vector  in  the  x,  y,  and  z  directions.  We  have 


=vi+vj+vk 


and  so 
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?  ?  ? 

V  -t  —  T  V 

x  y  z 


i  -*  j  2; 

t  v  i 


l  .  !j  'l 
\  -*-  *  / 


Since  {  v}2  is  xnovn,  Equation  (l4)  is  one  condition  for  the  determination  of  tie 
urknovns  v^,  v„,  v^.  We  need  tvo  more  conditions.  We  tfiil  assume  t  .at  the  vector 
v  lies  in  the"  o  lane  determined  by  'u  and  n  (see  Ti.-ore  4).  ibis  assumption  can 
be  stated 


/  -i  v  » 

(  u  x  n}.v  =  0 


(15) 


Cince 


a  -» 

u  x  n 


u  u  0 

*  y 

n  n  n 
x  y  z 


and  since  v  is  given  by  Equation  (13),  Equation  (15)  reduces  to 

v  n  v  -  v  11  u  ♦  v  (u  n  -  u  n  )  =  0 
xzy  yzx  z  xy  y  x 


(16) 


■deferring  to  Figure  4,  we  have 

u  .  v  =  I  u|  i"v|  cos(6  +  0) 


(17) 


and  since  u.v  =  u  v  +u  v  ,  we  have  from  Equation  (17) 

x  x  y  y 


vu  tv  U  =  |v|  cos  (0  +  0) 

xxyy 


(13) 


'Trip  onlv  unknowns  Ip  Koiiat.lons  (in)  ana  (lS)  are  v  ,  v  ,  v  :  and  (l6)  ar.a  (l5) 
w  -  x  y 

are  linear  in  v^,  v  ,  v^.  Hence,  it  is  easy  to  solve  (16)  and  (l3)  simultaneously 
for  v  and  v  in  terms  of  v  and  known  quantities.  The  results  are  as  follows: 

x  y  z 


Vz(Vx  "  W>  *  Vz  M  cos  (Q  +  p) 

l-f.  2 
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/,  ~\ 
u*/ 


V  5= 
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nzuv  |u|  !  v|  cos(0  •»  0)  -  uvz  (uynx  -  uxny) 


{  "A  '• 
\^W 


n  lul 
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and  Plane  detemined  by  u  and  n 


FIGURE  4. 


Substituting  Equations  (±9)  and  (20)  for  v„  and  v  in  Equation  ( l-'t ) ,  we  obtain 
an  equation  in  which  the  only  unknown  if>  v  .  This  equation  is  easily  solved  to 
give 

(21 


v  and  v  can  now  be  computed  by  substituting  this  solution  for  v  into  Eauaxio 

^  y  .  . 

(19)  and  (20).  Finally,  we  have  (see  Figure  3)! 
a  =  sin-1  (Vy/  |v|  ) 

=  tan"1  (v,/  vx) 

III.  ESTIMATION  OF  DANGER  AREAS 

The  theory  developed  in  Section  II  can  be  applied  to  estimate  danger  areas 
due  to  ricocheting  projectiles.  For  the  particular  projectile,  gun  and  line  of 
eire  under  consideration,  the  velocity  of  impact  ju]  and  angle  of  fall  are  known 
or  can  be  computed  for  any  range,  R.  Consider  a  particular  range,  il.  We 
suppose  the  projectile  strikes  an  inclined  plane  at  R  .  As  the  orientation 
of  the  inclined  plane  is  varied,  by  varying  the  values  of  n,  and  n  ,  the 
various  ricochet  ranges  are  determined  by  means  of  the  theory  of  Section  II  and 
the  rest  coordinates  after  ricochet  are  plotted  in  the  x-s  plane.  This  procedur 

is  repeated  for  other  pertinent  values  of  R,  say  Rg,  R^,  . R^  ,  and  for  each 

range  the  set  of  rest  coordinates  is  plotted.  An  envelope  of  the  family  of  set 
points  is  sketched  in  to  form  the  boundary  of  the  danger  area.  In  the  interest 
of  safety,  it  is  generally  assumed  that  the  projectile  fiys  just  sc  well 
after  ricochet  as  it  did  before  ricochet  —  i.e.  the  same  ballistics  are  used 
+.0  run  the  ricochet  trajectories  as  were  used  to  rua  the  original  trajectories. 
Cne  exception  to  this  is  that  windshields  are  assumed  to  break  off  after  impact 
on  hard  surfaces.  In  any  Ouse,  however,  it  is  assumed  that  the  flight  0/  the 
projectile  after  impact  is  stable. 


(22 

(23 
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An  example  of  a  danger  Rrea,  constructed  according  to  the  shove  procedure, 
is  given  in  Figure  s. 

IV.  VERTICAL  DANCER  AREA 

In  certain  cases.  a  training  area  may  be  located  in  the  vicinity  of  an  airfield 
and  training  exercises  may  be  in  progress  while  aircraft  are  taking  off  or  landing. 
It  is  of  interest  to  determine  the  veilb-a.1  area  dangerous-  *o  the  aircraft  because 

of  ricocheting  projectiles.  This  problem  4~  a  special  case  of  the  lateral  ricochet 
problem  already  described  except  for  the  following  specializations.  The  values 
of  the  ricochet  variables  'v  and  ' t>  are  used  to  compute  the  2enith  coordinates, 
instead  of  the  rest  coordinates,  ofter  ricochet.  As  before,  the  coordinates  are 
c^puted  for  each  location  end  orientation  of  the  ricochet  plane.  The  totality  of 
the  zenith  coordinates  fill  a  3  dimensional  volume  whose  envelope  defines  +hp 
boundary  of  the  vertical  danger  area.  The  boundary  is  now  a  suriace  instead  of  a 
line.  The  section  of  the  envelope  in  the  vertical  plane  which  contains  the  line 
of  fire  is  ordinarily  of  most  interest  since  this  section  includes  the  highest 
zenith  points. 


u  ) 

D.  J.  DUNN 


VJ.  G.  DOTSON,  JR. 
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FIG.  5 

LATEHAL  RICOCHET  'PEA  OF  THE  (05MM  M392 
APDS-T  PROJECTILE  OFF  A  SAND  OR  CLAY  SURFACE 

OA  =  LINE  OF  FIRE 
OABCO  »  POSSIBLE  RICOCHET  AREA 


NOTEi  THE  RICOCHET  AREA  SHOWN  LIES  TO 
THE  RIGHT  OF  THE  RIGHT  LIMrr  LINE  OF 
FIRE.  A  SIMILAR  AREA  LIES  TO  THE  LEF1 
OF  THE  LEFT  LIMIT  LINE  OF  FIRE. 


V4  REFERENCES 


*,  TJ-I 

■L  .  « 


pt<  TmiH  Surfdoss.  rrt.  Ssport  No*  Moy^V'  i oii^ 


2.  Hitchcock,  H.  P.  Effects  of  Ricochet  on  the  Motion  of  ‘*1  oj^ctiles.  BRL 
Report  No.  62^,  February  lQh-7, 


DISTRIBUTION  LIST 


of 

Organization 


20  Commander 

Defense  Documentation  Center 
ATTN:  TIFCR 
Cameron  Station 
Alexandria,  Virginia  22 Jilt 

2  Commanding  *** 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD-RP-B 
■Washington,  D.  C.  20315 

1  Commanding  Officer 
Frankfora  Arsenal 
Philadelphia,  Pennsylvania  19137 

1  Commanding  Officer 
Ficatinny  A.-  seuai 
Dover,  New  Jersey  O70O1 

1  Commanding  General 

U.  S.  Army  Missile  Command 
Redstone  Arsenal,  Alabama  35809 

1  Commanding  General 

U.S.  Army  Munitions  Command 
Dover,  New  Jersey  O78OI 

1  Commanding  General 

U.  S.  Army  Weapons  Command 
Rock  Island,  Illinois  61200 


Crojes  Organi  zation 


1  Commanding  Of l rcer 

U.  0.  Army  Artillery  Combat 
Developments  Agency 
Fort  Sill,  Oklahoma  7350J 

1  Commanding  Officer 

U.  S.  Army  Infantry  CD  Agency 
Fort  Benning,  Georgia  31905 

1  Commanding  Officer 

U.S.  Army  Armor  CD  Agency 
Fort  Kncx,  Kentucky  40120 

1  Commanding  Officer 

U,S.  Army  Aviazion  CD  Agency 
Fort  Rucker,  Alabama  56562 

1  Commanding  Officer 

U.S.  Army  Special  Warfare  CD  Agenc 
Fort  Belvoir,  Virginia  22060 

1  Commanding  General 

U.  S.  Army  Combat  Developments 
Center 

Fort  Ord,  California  95941 

2  Commanding  General 

U.  S.  Army  Continental  Army 
Command 

Fort  Monroe,  Virginia  23551 


1  Commanding  General 

U.  S.  Army  Electronics  Command 
Fort  Monmouth,  New  Jersey  077 ~ * 


1  President 

U.  S.  Army  Infantry  Board 
Fort  Benning,  Georgia  31905 


2  Commanding  General 

U.  S.  Army  Combat  Developments 
Command 
ATTN:  CDCMR 

CDCRE-C 

Fort  Belvoir,  Virginia  22060 


1  President 

U.  S.  Army  Armor  Board 
Fort.  Knox,  Kentucky  40120 

2  President 

U.  S.  Army  Artillery  Board 
Fort  Sill,  Oklahoma  '75503 


1  Commanding  General 

U.  S.  Army  Combined  Arms  Combat 
Developments  Agency 
Fort  Leavenworth,  Kansas  66027 


BISTRIPt/TT'JM  I.IST 


No.  of  No.  of 

Copies  Organization  Copies  0-  gas 'ration 


2  Commandant 

U.  S.  Army  Infantry  School 
Fort  Benning,  Georgia  3190 5 

2  Commandant 

U4  ft  A  y»py  Avwqw  CCV,0G1 

Fort  Knox,  Kentucky  +0120 
5  Commandant 

U.  S.  Army  Artillery  and  Guided 
Missile  School 
Fort  Sill,  Oklahoma  73503 

1  John  F.  Kennedy  Space  Center, 

MSA 

ATTN:  Library 

Cocoa  Beac+i  Florida  5295-1 

1  Commander 

U.  S.  Naval  Ordnance  Test  Station 

ATTN :  Library 

China  Lake,  California 

3  Chief,  Bureau  of  Naval  Weapons 
ATTN:  DL1-3 

Department  of  the  Navy 
Washington,  D.  C.  20360 

1  Commander 

U,  S.  Naval  Weapons  Laboratory 
Danigrea,  Virgi  a  22bi;8 

1  Director 

Marine  Corps  Landing  Force 
Development  Center 
Narine  Corps  School 
Quantico,  Virginia  2213b 


j  J  jLT'indfint 

Head quarters.  Marine  Corps 
Washing' on,  D.  C.  20380 

1‘  Australian  Group 

-/ o  A  r+a.-he 

2001  Connecticut  A. venue,  N.  W. 

Washington,  D .  o .  20Q0S 

10  The  Scientific  information  Officer 
Delence  Research  St.  " 

British  Embassy 

5100  Massachusetts  Avenue.  N.  W. 
'Washington,  D.  C.  20003 

Defence  Research  Member 
Canadian  Joint  Staff 
21+50  Massachusetts  Avenue,  N.  W, 
Washington,  D.  C.  2000b 

Aberd  on  Proving  Ground 

Chief,  TIB 

Air  Force  Liaison  Office 

w:  1UU  UViUO  sic.  isuli  ± 

Navy  Liaison  Office 
CBC  Liaison  Office 

D  k  PS  Branch  Library 

U.  S.  A.rmy  'Test  a  Evaluation 
Command 


18 


,**flW*#*«,H*  I'-B^^rfifs  v’WW'*-  mf|  i>WW^Wnri|‘»’'  W“‘ 


